Free Expansion of a Weakly- interacting Dipolar Fermi Gas 
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Abstract 

We theoretically investigate a polarized dipolar Fermi gas in free expansion. The inter-particle 
dipolar interaction deforms phase-space distribution in trap and also in the expansion. We exactly 
predict the minimal quadrupole deformation in the expansion for the high-temperature Maxwell- 
Boltzmann and zero-temperature Thomas-Fermi gases in the Hartree-Fock and Landau- Vlasov 
approaches. In conclusion, we provide a proper approach to develop the time-of-flight method for 
the weakly-interacting dipolar Fermi gas and also reveal a scaling law associated with the Liouville's 
theorem in the long-time behaviors of the both gases. 
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Recent development of trapping and manipulating techniques of atoms enables to realize 
ultracold dipolar gases of atoms and molecules in several years. In fact, dipolar Bose gases 
are realized in 52 Cr atoms with magnetic dipole moments and dipolar Fermi gases are 



also realized in heteronuclear 40 K- 87 Rb molecules with electric dipole moments [^-l6(. The 
dipolar gases offer a great infrastructure for quantum many-body physics as new types of 
actual quantum systems. 

The dipolar gases have a tensor- type dipolar interaction 7|, which depends on not only 
inter-particle relative position but also directions of the dipole moments. When we suppose 
the dipolar fermions perfectly polarized in a strong external electric field along the z-axis of 
the Cartesian axes as in the actual experiments |4t6(, the interaction becomes 

1 - 3r 2 Jr 2 



v(r) = G a 



(1) 



with r = {r x ,r y ,r z }, r = |r|, and a positive coupling constant Gd oc d 2 given later with 
some scaling units, where the electric dipolar moment d is experimentally manipulated in 
d ~ 10 _11 ~°eaB by choosing molecular internal states and electric field strength 4|-[6|]. 

As an important feature, the dipolar Fermi gas is purely dominated by the long-range 
tensor force because of vanishment of short-range interactions owing to the Pauli blocking 
effect. The axisymmetric interaction in eq. flTJ produces anisotropic quantum correlations 
and density deformation and also makes instability in strong interaction regime js i^j]. Fur- 
thermore, differently from the short-range interactions, the long-range dipolar interaction 
deforms phase-space distribution, which indicates important information for quantum many- 
body effects in the dipolar gas 

Free expansion physics is quite important for the dipolar gas. That is because the ex- 
pansion is often used to observe the momentum distribution of the trapped gas in the 
time-of-flight (TOF) method; and the anisotropic momentum deformation directly reflects 
the phase-space deformation and interaction effects in the dipolar gas. In fact, the dipolar 
Bose gases are firstly observed in the expansion lM3(. The dipolar interaction may change 
the momentum distribution not only in trap but also in the expansion, then we need some 
corrections to apply the TOF method to obtain the initial momentum distribution differ- 
ently from the ballistic expansion with a conserved momentum distribution. To obtain the 
corrections, some people develop the quadrupole scaling method [lo|, Ijj]; however no work 
has been done to demonstrate exact time-evolution of the expanding Fermi gas and confirm 
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validity of the scaling method in the expansion. The expansion problem is one of the most 
important current issues in theoretical physics on the dipolar gas. 

The aim of this work is to reveal the expansion problem, i.e. relationship between the ini- 
tial and final momentum deformations, in weak interaction regime. Here we exactly evaluate 
time-evolution of the minimal quadrupole deformation of the expanding gas by developing 
another scaling ansatz, which reproduces the Liouville's theorem and perturbation theory 
from the ballistic expansion. This work must provide a proper correction for the TOF 
method to detect the quantum many-body effects on the dipolar gas. 

Let us define the Hamiltonian of the dipolar gas as 

1. 



z/>(r) + H, (2) 



H = J dr ^(r 
with the fermion field operator ip and interaction part 

Hi= [ dr [ dr'i)\r)i) ] (r')v(r - r')i/;(r')ip(r), (3) 



where v is the dipolar interaction in eq. (CQ). The trap term is given as v^{r) = Q(—t)Vr(r) 
with the step function O(x) and trap potential Vr(r) before t — 0, when we leave the gas. 
To simplify the notations, we choose units for the reduced Planck constant and mass as 
H — 1 and m — 1 and omit the time parameter t, e.g. ip(r,t) = ip(r). 

To treat Hi in eq. ([2]), we apply the time-dependent Hartree-Fock approximation (TD- 
HFA) 



131 ]. in which the two-body interaction Hi is rewritten into one-body interactions with 
self-consistent mean-fields, the direct (Hartree) and exchange (Fock) terms, reflecting the 
many-body effects as mean-values. In TDHFA, the field operator satisfies an equation of 
motion, 

dijj(r) 
i- 



\v%^{r) + J dr'U(r,r')i;(r'), (4) 

consistent non-local se 

a 



dt 

with the self-consistent non-local self-energy U(r, r') given later in eq. (|7J). Then the Wigner 
function 

/(r,p) = /rf S e-^ s n(r-|,r + |) (5) 
with the density matrix n(r, r') = (^(r')ijj(r)), obeys 



with the Wigner transformation of U(r, r'), 

u(r,p) = Jdse- i P- a u(r + ^,r-~ 

= v T (r) +u u (r) +u F (r,p), (7) 

where and V" act on / and u, respectively. Here we introduce the Hartree part 

ii H (r) = / dr'v(r — r')n(r') (8) 



with the number density n(r) = n(r, r) and Fock part 

dq 

(27) 

with the Fourier transformation of v(r) in eq. (CQ) 



u*(r,p) = -J ^V{p-q)f{r,q) (9) 



V(p) =Jdr e-V*v{r) = -^G d ^-^, (10) 



where p = \p\. We explicitly write H (= 1) only in eq. (JB} for the later explanation. 



To treat t 
approach 



re quantum dynamics in eq. flSJ), we use the semi-classical Landau- Vlasov (LV) 
15| ; i.e. expanding the trigonometric function in eq. ([6]) up to the first order of h. 
Then we obtain the LV equation, 

^ = - (p + v pU ) ■ (V r f) + (Vr«) • (V P /) . (11) 

As a conclusion of the formulation, we study the expansion by solving the LV equation, 
eq. (ITT]), with some initial conditions, f(r,p;t = 0) = f (r,p). 

Here we emphasize that a feature of the finite-range interaction (e.g. eq. ([TJ) must be in 
contribution of the Fock part, i.e. the momentum dependence of u as shown in eq. (JTJ), so 
that VpM in eq. ( II lj) may give notable contribution in the dipolar gas and should not be 
omitted. 

In this work, we focus on the interaction effects on the time-dependent momentum distri- 
bution, p(p) = J drf(r,p), which corresponds to the density distribution at long-time limit 

as 

"M=(2^(l) +0(rl >- (12) 

That is because we can pragmatically neglect the interaction after a cut-off time t c (— >■ oo in 
this theoretical work) [16[ , when the Wigner function in eq. (II lj) is regarded as the ballistic 
solution, 

f(r,p) = f (Z,p), (13) 



with the Galilei transformation, r — > £ = r — pt; then the density distribution at t > t c 
becomes 

dp 



n(r) = J -^f c (r-p(t-t c ),p) 



r/c^ + ptcyl+Or 1 ) (14) 



(2?rt) 

with the Wigner function f c (r,p) at t = t c , where eq. ( I14p agrees with eq. (I12p owing to 
the conservation of the momentum distribution in the ballistic expansion at t > t c . Thus 
the interaction effects appear in the time-evolution before t c , and one can obtain the final 
momentum distribution from the density distribution at long-time limit as shown in eq. (fl2j) . 

Before calculation of the expansion dynamics, we firstly consider the ballistic expansion; 
i.e. the interaction is neglected only in the expansion process except in the initial states. In 
the ballistic expansion, the LV equation, eq. (fTTj) . becomes 

f = -p-Vr/. (15) 
Its solution exactly agrees with the ballistic solution in eq. f fl3|) . When we assume the initial 



Wigner function as some function of a parameter 13] 



*{r,p)= £ + (16) 

l=x,y,z 

the dynamics is determined only by a time-dependent parameter z = z(£,p), which becomes 
z = ufpft 2 — 2uifrip t t + 0(t°) at the long-time limit. 

In the previous approach [12J, the dynamical parameter is approximated as z ~ 
z(r',p') with r[ = b^ri and p[ = bipi — b^i, where bi obeys bi = Dfbf 3 with Di = 
y/(p})(t = 0)/(rf)(t = 0) and b t (t = 0) = 1 and becomes b t = Dit + l/(2Dit) + 0(t -3 ) at the 
long-time limit, after all. Here the cap-dot represents the Newton's notation for time deriva- 
tive. Owing to the quadratic Wigner function, D\ = uii/ji, and z ~ ufpft 2 — 2ufripit + 0{t°) 
in the previous ansatz. This result agrees with that of the exact solution. However this 
agreement is a specific property of the quadratic Wigner function. Our ansatz in this work 
generally satisfies the ballistic limit beyond the quadratic case and gives the exact solution 
for the minimal variation from the ballistic expansion as described below. 

We now consider the quadrupole deformation of the momentum distribution in the in- 
teracting non-ballistic expansion. The quadrupole deformation is the most important de- 
formation owing to the axisymmetric interaction in eq. (00) and reflects not only the initial 



deformation in trap but also the additional deformation in the expansion. To describe 
them, we introduce a time- dependent index parameter: X p = ln(7^/7o) with T z = (p 2 z ) and 
7o = {(fix) fay) (jFz)) ■> then our interest is just focused on relationship between the initial 
conditions and X p (t — > oo). Eq. f lTTj) leads 

^ " ^ f^| ^ Mp), (17) 



dt J J (2tt)3 T z d£ z 

where we define f(£,p) = f($ + pt,p) and use the partial integral. According to the 
definition in eq. (J7J), u(£ + pt,p) in eq. (TTTj) can be written as 

u(£ + pt,p) = | dry | t^K* - V, (P ~ <z)*)/(*7, q), (18) 
where we introduce 

v(x, at) = v(x + at) - V(a)5(x + at) (19) 

with q = a + p and rj = x + By substituting eq. (fl8"|) into eq. ( TIT)) and using exchange 
symmetry on {£,p} -H- {r?,^}, we obtain 

rfAp - ^ [fo^( Xi g i ) dF ( x > a! ) (20) 



dt J J (2tt) 3 T z y ' ' dx 
with 

F(*,a) = JdtJj0-m,p)f(Z + x,p + a). (21) 

We here assume small variation from the ballistic solution in the phase-space distribution 
and introduce an ansatz for the Wigner function as 

Kt,P) = Mlp) (22) 

with 

Px, y = e x/A (p x , y + j£ x , y ), p z = e~ x/2 (p z - ~£ z ), 

£,x,y = 6 ^ ^>x,yi £,z = & ^ £,zj (23) 

where A indicates the additional deformation in the expansion, X p ~ Ao + A, with X(t = 0) = 
0. At the ballistic limit, A = 0, and it always reproduces the ballistic solution. Furthermore 
this ansatz exactly satisfies the Liouville's theorem corresponding to the current equation in 
the microscopic theory Ji], Q . 
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In this Letter, we consider two general initial conditions: the high-temperature Maxwell- 
Boltzmann (MB) and zero-temperature Thomas-Fermi (TF) gases in cylindrical harmonic 
oscillator traps with the small phase-space quadrupole deformation owing to the weak in- 
teraction. The initial Wigner function becomes 

f (r,p) = e-( RM+p < +p ^ 2 (24) 

for the MB gas and 



/o(r,p) = e 

for the TF gas with 



_ R c + R z + P c + P z 



(25) 



R c = e^- x ^ 2 (r 2 x + r 2 y ), R z = e^-^r 2 , 

P^e^ipl+pl), P z = e~^vl (26) 

where Ao = \ p (t = 0) indicates the initial momentum deformation {2(| , and Ao indicates 
sum of the initial density and momentum deformations. In addition, we also determine 
Gd in eq. ([T]) as G d = (k B T/(hu))-^ 2 N(a B /a ho )(m/m ( Md/(ea B )) 2 for the MB gas and 
Gd = (6iV) 1 / 6 (aB/aho)(^/'^e)(^/(eaB)) 2 for the TF gas [21| with the Boltzmann constant 
k B , temperature T, particle number N, Bohr radius Ob, electron mass m e , elementary electric 
charge e, trap frequency u = (tOxCOyUz) 1 ^ 3 (= 1 in this Letter), and oscillator length a^o = 
%/{muS). The above formulations reproduce general quadratic descriptions of the MB and 



TF gases with rescaling 



22| 



Here we should comment that the interaction effect on A in eq. (|26|) is always cancelled 
out in the quadrupole deformation owing to conservation of phase-space volume in the 
Liouville's theorem. Thus A corresponds to the density deformation of a non-interacting 
gas in the initial trap potential, i.e. Ao = In (u I /u 2 ) 4 '' 3 = In (tOy/tOz) 4 ^ 3 . 

According to eqs. fT2"2l . fT24|) . and f )25|) . F(x,a) in eq. fl2T|) becomes 

F(x, a) = V* 2 / 4 (27) 

for the MB gas and 



for the TF gas, where 



s = JE Aj+A M + £ X j + X z (29) 

with 



3=X,V J=X,r 



Substituting eq. (|27|) or ( 128|) into eq. (120]) and expanding up to the first orders of A , A p , 
and A, we obtain a time-evolution equation for the minimal deformation: 

oo oo oo 

A = G d E E ^ukK Ap A^ 

j=0 J=0 i^=0 

— C d (r 000 + TioqAo + r 010 A p + TooiA) (31) 

with the time-dependent coefficients T/j^(t) depending on the initial conditions. 
After the calculation, we obtain r 00 o = r io = r 001 = and 

A ^ G d A r 100 (32) 

in the both gases. It reveals that the additional deformation A in the expansion depends 
only on a conserved quantity G^Ao = S and exhibits a scaling behavior, 

X(t ->• oo) ~ 7i 00 G' d Ao = TiocS 1 , (33) 

at long-time limit with 7100 = Jo°° ^ -^100- The scaling behavior is due to cancellation of 
the Hartree and Fock parts determined by the time-developing phase-space distribution 
with the quadrupole deformation. According to the Liouville's theorem, the incompressible 
quadrupole deformation in the momentum space gives strong constraint for the spacial 
motion as shown in the description of the deformation ansatz, and this constraint provides 
a source of the cancellation and scaling behavior. 

Here it should be noted that the minimal deformation in eq. fl32l) exactly agrees with that 
in the first order perturbation theory, which can be obtained by replacing F(x, a) defined in 
eq. (T2"U1) with that of the ballistic solution, i.e. substituting A = and A = in eqs. (1271) and 
(1281) . This agreement is obviously due to the cancellation in the non-perturbative relation 
in eq. ( 131]) and, after all, indicates that the minimal deformation must be determined only 
by the first order perturbation according to the scaling behavior. 
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In addition, we obtain 

ri0 ° = (f +1)5/2 ( 34 ) 

with 7 100 = -3/(70^) « -0.0242 for the MB gas and 7100 ~ -0.00444 for the TF gas. 
Thus the quantum effect, i.e. difference between the MB and TF gases, is in the values 
of 7100 and definitions of Gd- According to the temperature dependence in the MB gas 
(Gd oc T -5 / 2 ), the deformation vanishes at the high-temperature limit; in other words, the 
interaction effect principally appears as the quantum effect. 

As a result, according to eq. fl32|) and Gd > 0, the interaction slightly reduces (increases) 
\p when A > (< 0). It is due to the angular dependence of the dipolar interaction in 
eq. dl]). When A = 0, corresponding to the spherical trap, the results reproduce those of 
the ballistic solution, A = 0, owing to vanishment of the interaction effect by the angular 
integration. 

Note that S in eq. (133]) can be given by the experimental setup, and A p (oo) (~ Ao + A(oo)) 
can be measured from the density distribution at long-time limit according to eq. ({12]) . Thus 
eq. fl33|) provides a theoretical approach to obtain the initial momentum deformation A 
from the experimental measurement as a proper correction for the TOF method. Because 
of smallness of 7100 for the both gases, the correction is negligible when S (or the trap 
anisotropy) is small; otherwise the correction must have important contribution to detect 
the interaction effects. 

Finally we should comment on difference between our and the previous approaches 



10. 
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to treat the quadrupole deformation in the expansion. These are theoretically same except 
for choice of the phase-space frames: the expanding frame in this work and rest 

frame {r,p} in the previous works. Thus any difference in the results must be due to 
the frame transformation. The previous ansatz is originally developed to the expanding 
hydrodynamical Bose gases, and then reproduces agreeable results with the experimental 
measurements 23], 24]. On the other hand, our ansatz is made to reproduce the exact results 
for the minimal deformation from the ballistic solution in the weak interaction regime and 
give explicit description of the variation. The both approaches give same results at the 
ballistic limit for the quadratic Wigner function as discussed above. Thus we find out that 
the previous ansatz also reproduce our exact results up to the first order perturbation from 
the ballistic limit in the quadratic case. That is because the exact results depend only on the 
first order perturbation according to the scaling behavior, and the first order perturbation 



is determined by the ballistic solution according to the perturbation theory. This nontrivial 
agreement may be a significant property of the quadratic deformed gas associated with the 
scaling law. 

In conclusion, we reveal that the additional phase-space deformation in the expansion is 
negligible when S is small, as shown in eq. f !33p ; as S increases, the deformation linearly 
grows and exhibits the scaling law associated with the Liouville's theorem in the long- 
time behaviors of the MB and TF gases. After all, our results give the exact prediction 
of the minimal quadrupole deformation in TDHFA and the LV approach, and provide a 
proper correction for the TOF method as a fundamental infrastructure for the dipolar gas 
physics. At the end, it should be noticed that the scaling law may not be valid for the large 
deformation by the strong interaction; then one should directly solve the time-evolution in 
the LV equation beyond the quadrupole ansatz to include the other multi-pole contribution. 
The expansion problem in the strong interaction regime should be studied in another paper. 
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